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The phase dynamics of a long Josephson junction with spatially inhomogeneously distributed 
bias current is considered for the case of a dense soliton chain (regime of the Flux Flow oscillator). 
To derive the analytical solution of the corresponding sine-Gordon equation the Poincare method 
has been used. In the range of the validity of the theory good coincidence between analytically 
derived and numerically computed current-voltage characteristics have been demonstrated for the 
simplest example of unitstep function distribution of bias current (unbiased tail). It is shown, that 
for the considered example of bias current distribution, there is an optimal length of unbiased tail 
that maximizes the amplitude of the main harmonic and minimizes the dynamical resistance (thus 
leading to reduction of a linewidth). 
PACS number(s): 74.50. + r 



I. INTRODUCTION 

Long Josephson oscillators operating in the flux-flow regimeS are-, presently considered as possible 
devices for applications in superconducting millimeter-wave electronicsa,B. In comparison to single fluxon 
oscillators they have higher output power, wider bandwidt 
wider linewidth of the emitted radiation from the junctional 
junctions has been intensively studied both numerically and analytically. As classified inQ, the studies can 
be divided into two categories. In the first category one can find sine-Gordon solitons created by a small 
magnetic field confined to the boundaries of the junction and propagating along the junction under the 
influence of the bias current. In this case one can extract detailed information about the physics of the 
problem from the McLaughlin and Scott theoryu as the aim of that theory was to investigate the steady- 
state motion of fluxons (and other solutions of the sine-Gordon equation) under small perturbations. In 
the second category one can find dense soliton chains distributed along one dimension of the junctions 
driven by dc and (or) ac currents. The physics of this category cannot be treated by the McLaughlin-Scott 
perturbation theory for solitons as the perturbations (bias current and magnetic field) are no longer small 
compared -ta-the sine-Gordon terms. However, namely this regime is interesting for practical flux flow 
oscillatorsEl,0,ElJa. The latter regime is called the "flux-flow" regime and is characterized by excitations 
which travel on top of a fast rotating background so that the effective nonlinearity in the system is 
drastically reduced due to fulfilling the following conditions: 770/0; 2> 1 and (or) r » 1 (where 7]$, a 
and r arc dimcnsionlcss total dc bias current through the junction, the damping and the magnetic field, 
respectively). It should be noted, that the "flux-flow" regime has nothing to do with motion of well- 
distinguishable flux-quanta: the soliton chain is so dense, that one should speak about transmission of 
quasi-linear waves in a long Josephson junction in this case. n 

For the flux-flow-cegime two main approaches are known. One is the approach by KulikB (that has also 
been used, e.g., ir£3,tj), that is based on linear mode theory and perturbative analysis around rotating 
background (<j) = 4>o + ip, ip <C 1). Another approach has been suggested in Ref.ta and is also based 
on known form of a solution around which one can make a perturbative expansion, but in a different 
manner: the unidirectional fluxon train is accompanied by two plasma waves, that allows to satisfy 
boundary conditions. Using either of the approaches one can compute the current- voltage characteristic 
of FFO in the second order approximation for spatially homogeneous dc and ac driving. However, if 
the bias current has inhomogeneous spatial distribution, it is not clear how to derive the corresponding 
solutions, since the approaches are based on the "anzats", some k-apwn form of the initial solution, around 
which the perturbative analysis has to be done. It is known fromll3 that the use of inhomogeneous driving 
(unbiased tail) may lead to decrease of a dynamical resistance and therefore to reduction of a linewidth. 
Also recently importance of accounting-af magnetic field fluctuations for the analysis and design of FFO 
has been experimentally demonstrated!!^. Further theoreticaLJ and experimentalli3 investigations have 
shown that not only usual dynamical resistance Rd = dV/dl, but also the dynamical resistance with 
respect to the magnetic field = dV/dH (that may also be attributed to the dynamical resistance of 
control line R^ L = dVjdlch) is important for calculation of the linewidth of FFO. Therefore in FFO 
designs one should care about minimization of both above mentioned dynamical resistances and the 
present theoretical state of the problem should be reconsidered in order to improve characteristics of 
practical flux flow oscillators since nobody payed attention before to the value of R^ . 



The aim of this paper is to present the approach that allows to systematically study the regime of 
large magnetic fields and bias currents and to describe the dynamics of the Flux Flow oscillator with 
spatially inhomogeneous driving (and, as particular example, with unbiased tail). This approach allows 
to perform the detailed analysis of FFO and to find an optimal profile of the bias current distribution 
and an optimal length of unbiased tail in order to maximize the output power of FFO and to reduce the 
lincwidth. 



II. BASIC EQUATIONS AND THE POINCARE METHOD 

The electrodynamics of a long Josephson junction in the presence of magnetic field is described by the 
perturbed sine-Gordon equation 



<9 2 <90 d 2 
dt 2 dt dx 



a S7 ~ 7^2 = ^( x ) ~ sin (0) C 1 ) 



subject to the boundary conditions 



d<t>(0,t) = d<t>(L,t) = r 
dx dx 



In this equation space and time have been normalized to the Josephson penetration length Aj and to 
the inverse plasma frequency w~ , respectively, a is the loss parameter, rj(x) is the normalized dc bias 
current density and T is the normalized magnetic field. In accordance with RSJ modeOjla one takes the 
loss parameter a = — , where u> p = ^/2eI c /hC, u> c = 2eI c RN /%, C is the capacitance, Rn is the normal 

state resistance (i?jv = V/I qp , V being voltage and I qp - the quasiparticle component of the current), 

i i 

I c is the critical current, rj(x) = J(x)/J c (I = J J(x)dx, I c = J J c (x)dx, I is the bias current), I is the 

o o 

dimensional length of the junction, L — l/Xj. 

As it has been mentioned in the introduction, the flux-flow regime is characterized by the fulfilling of 
the following conditions: 770/a 3> 1 and (or) r» 1, where 770 = J rj(x)dx/L is the dimensionless total 
dc bias current. nr~l n 

Instead of linear mode theory_and perturbative analysis around rotating background@J13,EJ, one can use 
more general Poincare methodll3: obtain the solution as the series with respect to the naturally arising 

in the flux-flow regime small parameter e = (^-^j *C 1. Let us change variables in Eq. ([!]), r = 2e.£, 

z='^x: 

a 

d 2 ^ d<t> d 2 4> (3 . 

d^ + ^-d^ = ^ {z) - esm{ ^ (3) 

where (3 — a 2 /r/ . 

The steady-state solution of this equation may be found in the form: 4>(t) — (/>o( r )+e0i(r)+e 2 02(T) + . . . 
(|0o (t) I ^> €\4>\{t)\ e 2 1 02 (t - ) I ^ ■•■)• Substituting this into Eq. (||) one can find the zero order 
equation: 

dr z dr dz z rja 

It is easy to see, that the steady-state solution of this equation is: 4>q(t) — r + jz — g(z) — ^t + Tx — g(x), 
7 = aT/rjQ, where: 

g(x) = ^ —— / r)(x) cos k m xdx, (5) 

m=l Lkm Jo 

k m = TTm/L, so in the 0-order approximation the current-voltage characteristic is given by the ohmic 
line: v(r]o) = d<f>/dt = flj(r/o) = r/o/a (due to Josephson relation the voltage is proportional to the 
Josephson frequency j, here means averaging in time). To get higher order equations let us decompose 
sin(0o(r, z) + e0i(r, z) + e 2 02(r, z) + . . .) into Taylor expansion. From the structure of the considered 
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linear recurrent equations it is clear that the steady-state solution </> n (r, z) may be presented in the form: 
4>u{t, z) = uj n T + <pn P (T, z), where 0„ p (t, z) is periodic nongrowing component. 

Collecting together all linearly growing components w„t one can get: sin(0(r, z)) = sin({wo r + ew]T + 
e 2 LL> 2 T+- ■ .+jz — g(z)} + e4>i p (T) + e 2 (l)2p(T) + . . .). Now one can linearize sin(0) as: sin(</>) w sin(ujjT+jz — 
g(z))+e(j)ip(T, z) cos(aj,/r+7z-g(z))+e 2 02p(r, z) cos(ujjT+jz-g(z)) + . . ., where uij = ui +€ui 1 +€ 2 ui2 + - ■ ■ 
is the oscillation frequency (u>q — 1), and w 2 ,..., w n ,..., 0i p (t, z), 4>2 P (t, z), </>„ p (t, z), n . r -are unknown 
functions that is required to obtain. One can consider the solution up to the 4-th order (inB-liil the solution 
up to the 2-nd order has been derived, but in principle it can be done up to any order, all equations may 
be solved recursively), and then the following equations for 4>\{t, z) - 04 (r, z) may be written: 

d 2 4>i , a d<f>i <9 2 </>i . 

'"~a = -sm{(jjT + yz- g(z)), (6) 



dr 2 8t dz 2 



d 2 4>2 




d 2 h 


dr 2 


dz 2 


d 2 fo 




d 2 fo 


dr 2 


dz 2 


d 2 <j> 4 




d 2 <f> 4 


dr 2 


dz 2 



-0i p (r) cos(cjjt + 7Z - g(z)), (7) 
-02 P (t)cos(cjjt + 7Z - g(z)), (8) 
-<f>3 P {T)cos(ujjT + -fz- g(zj). (9) 



The boundary conditions (p|) are taken into account in the solution of the zero order equation and Eqs. 

0. The solutions of linear Eqs. (p|)-(p|) 



@"@ should be solved for zero boundary conditions: |^ 



=0,i 



may easily be derived recursively substituting the solution in the form: </>„(r, z) = ^2 C nm {r) cosk m z. 

m=0 

The computer simulations of Eq. ([!]) have been performed on the basis of an implicit difference scheme 
(similar to one, presented in the Appendix oftij) with adaptively varying time step and calculation time. 



III. THE 4-TH ORDER APPROXIMATION FOR THE HOMOGENEOUS CASE 

Since the peculiarities of the derivation of the 2-nd order approximation (solution of equations (||),(0)) 
for homogeneous case have been considered in detail inla, the only solution given in original notations 
of Eqs. (0),(§) is presented below. The current-voltage characteristic ^(770) = Qj(r]o) = dcj>/dt may be 
obtained in the form: Qj^r/o) = r/o/a + ^(vo) + ^4(^0), where ^2(^0) an d ^4(770) are the second and 
the fourth order corrections of frequency (voltage), respectively. As it has been demonstrated inE_3, one 
has to solve transcendental equations for Q2 and f^4 in order to derive the voltage Oj as function of the 
bias current 770- However, in homogeneous case one can express the bias current ijo as function of £!j: 
770 (^2,/) = ol{Qj — Q,2(Qj) — VL/±(Tlj) — . . .} that gives an analytical expression for rjoijlj) and allows to 
avoid the solution of the transcendental equation. 

For the homogeneous case the 2-nd order correction as function of Q,j may be presented in the form: 

71=0 

_ rX(l - cos(rL) cos(Trn)) _ TL sin(TL) cos(Trn) 

Sn ~ {TL) 2 - (nn) 2 ' Cn ~ (TL) 2 - (7m) 2 

Correspondingly, the 4-th order correction may be written as: 

1 00 

^4 = — — ^ [Asulci - B 3 iilsi] , 

i=0 



(an,j) 2 + [k 2 



I- 



Cn\ 



(10) 



(11) 

(12) 



A 31i 



ail j [MimJs 

0(),i m=0 



B22m,Jcmi] ~ (k 2 — fl 2 ) \A<X2m,Jcmi ~ B22mJ Smi] 

m=0 



(aOj) 2 + (k 2 - n 2 ) 2 



Ban 



OtQ,J \Al1mJcmi — B22mJSmi] + {kf — fl j) [^22mJSmi + BllrnJcr, 

«0,i m=0 rn=0 



2 \2 



2 - So- 



A22m — 



_ 2 - S . 

t>22m — ~ — 



2aOj J2 [BlnJcnm ~ Ai n Jsnm) — (k m — 4f2 2 ) [AlnJcnrn + B\ n J Sum] 

n=0 n=0 



2aflj [AlnJCnm + Bi n Jsnm] + {k m — 4fij) [BlnJcnm — Ai n J$ 
n=0 71=0 



Ai n = (2 - 5 _) 



-Bin = — (2 — So,n) 



(2afij) 2 + (k 2 m - 4f>2)2 

anjicn - {kj - n 2 j)i Sn 
a nji Sn + (k\ - n 2 j)i Cn 



(anj)* + (kl ~ 02 )2 



1 J rLsinrLcos7r(m + n) ri sin TL cos 7r(m — n) 
Jcnm = ^ | (rL)2_ [^(m + n)] 2 + (r£) 2 - [7r(m - n)] 2 

1 J ri(l — cosTL cos7r(m + nj) TL(1 — cosrLcos7r(m — n)) 
Snm = 2 \ {YLf - [7r(m + n)] 2 + (TL) 2 - [n{m - n)} 2 

As one can guess, in order to obtain Jcmi and Jsmi, one has to interchange indexes in Jcnm and Js„ m , 
respectively; functions 7g„ and Ic n are given by (|ll|). 

The current-voltage characteristic ?/o(^j) = a {^J — ^2(^.7) — ^4(^7)} is presented in Fig. 1 for 
a = 0.1, L = 5, r = 3. It follows from the analysis that if the deviation from the ohmic line is located 
in the area where r)$/a > 2, it is usually enough to use the 2-nd order approximation only. It is seen, 
that the account of the fourth order correction gives more precise description of the height of Fiske steps. 
For description of Fiske or Eck steps that occur below the boundary t]q /a > 2 the 4-th order correction 
becomes of importance and significantly improves the approximation, see the inset of Fig. 1, where it 
is seen that the 4-th order approximation describes the step missed in the 2-nd order approximation 
(around f2j « 0.95). It should be noted that for larger magnetic fields the second order approximation 
coincides well with the results of computer simulation for rather small damping (e.g. for a — 0.04 — 0.01, 
see Fig. 2 for a = 0.04, T = 5, L = 5), but it is not easy to catch all Fiske steps in computer simulations 
in this case. 



IV. THE 2-ND ORDER APPROXIMATION FOR THE INHOMOGENEOUS CASE 

For inhomogeneous case, unfortunately, one can not avoid the solution of the transcendental equation 
and for simplicity let us consider the inhomogeneous case in the second order approximation only. For 
the second order correction Q2 of the current-voltage characteristic in the inhomogeneous case one can 
get the following transcendental equation: 



n=0 



((i 1 o+an2)) 2 + [k 2 ~( m /a + {l2) 



212 



(13) 



Isn = -j- J sin(Tx - g 



1 



L 



(x)) cos(k n x)dx, Icn = £ / cos(Tx — g(x)) cos(k n x)dx. (14) 







4 



It should be noted, that both in the homogeneous and inhomogeneous cases the equation for 0,2 is given 
by the same formula (|l^) with the only difference that in homogeneous case g(x) = in (|l4|) and functions 



Is n and Icn ma y be evaluated analytically (11) 



tail): rj(x) — rj s a(x — xq), rj s = ?° . In this case (k n = nn/L 



As an example of spatial distribution of the bias current let us consider the unit-step function (unbiased 

L—Xq ' 

2r) s _ i] s J (L - x Q )(x (2L - x ) - 3x 2 ), x < x , 



v ' ^ imkl 6L 1 x (2L z ~ 6xL + 3x z - xg), a; > a; v ; 

n— 1 v 

The transcendental equation (|l3|) with 5(0;) given by (|l^) may easily be solved and the voltage-current 
characteristic of FFO f2 ,7(770) = r/Q/a + ^iflo) niay be found. In Fig. 3 the results of computer simulation 
of Eq. (0) are presented for a = 0.3, L = 5, T = 5, —5, xq = 2.5 (long unbiased tail) and xq = (no 
unbiased tail). One can see, that in the case where the driving is homogeneous (xq = 0) the current- 
voltage characteristic is absolutely antisymmetric with respect to the bias current and the sign of magnetic 
field r does not play any role. Contrary, in the case with unbiased tail (xq — 2.5) the current-voltage 
characteristic becomes asymmetric: if the bias is applied at the end, where vortices exit the junction, the 
radiation is amplified, if the bias is applied at the end, where vortices enter the junction, the radiation 
is suppressed. The comparison between analytically derived IVC on the basis of Eq. (|l^) and results of 
computer simulation is presented in the inset of Fig. 3 for a = 0.3, L = 5, T = —5, x$ = 2.5. 

As it is seen from Fig. 4, not only the location of unbiased tail, but also its length plays an important 
role: the amplitude of the main harmonic has maximum as function of the length of unbiased tail (for the 
considered parameters the maximum is located around xq = 2.5; it follows from Eq. ([?]) that the second 
order correction of the IVC VI2 is proportional to the amplitude of the first harmonic <f>i(x,t)). In the 
inset of Fig. 4 the results of computer simulation of IVC minus ohmic part is given for different values 
of so and it is also seen that the main resonance is higher for xq — 2.5. It should be noted, however, 
that an optimal length of unbiased tail will be different for different parameters of the long Josephson 
junction (such as length, damping and magnetic field) and special investigations should be performed for 
optimization of designs of practical FFOs. 

It is clearly seen from Fig. 4 and the inset of Fig. 3, that the second order approximation of IVC in 
the inhomogeneous case gives rather good coincidence with the results of computer simulation, but the 
discrepancy increases with the increase of the length of unbiased tail, which is due to the fact that the 
deviation of IVC from ohmic line occurs for larger values of a/770 and the small parameter e = (a/770) 2 
becomes larger leading to worse applicability of the approximation. 

Let us perform the analysis of voltage-field characteristic of FFO to study the dependence of R*[ = 
dVj dH (or = dvj dT in dimensionless notation) on the length of unbiased tail. From the plots of fl j(T) , 
presented in Fig. 5, we do not see so significant dependence of from the length of the unbiased tail 
as it was for Rd- For smaller bias current the curve Oj(r) consists from several branches with positive 
derivative. However, for larger currents 770 > 1 it is intriguing to see that there is a pjige of parameters, 
where the derivative (Rf) takes negative values. As it has been demonstrated ir£3, if fluctuations of 
bias current rip{t) and magnetic field Tp(t) are correlated (such that Fp(t) — arjplt), a is a numeric 
coefficient), the linewidth Af of FFO is proportional to (Rd + oR^tf (where R d = dV/dl, Rf = dV/dH) 
and if Rd < oR^ one can see the plato on the plot of Af(Rd) (seetil), that is, with further decrease of Rd 
the linewidth does not change. It is obvious that if Rd > and crRd < (meaning that the noises of bias 
current and magnetic field are anticorrelated) the linewidth may be significantly decreased. Certainly, 
in real situations even if \Rd\ = \&Rd\ ( or i alternatively, if Rd and crR^ are both positive, but very 
small), the linewidth willjifiver be zero since other noise sources (e.g., technical fluctuations) that were 
not taken into account inEEl would become of importance and the saturation of Af(Rd, R¥) should be 
observed. It follows from estimations that FFO operating in the considered range of parameters (small 
length and large damping) will have rather small power that may be not enough for practical applications. 
Nevertheless, since there is still a lack of understanding about nature of magnetic field fluctuations, it 
could be interesting to experimentally study this case for ctR^ < to check if bias current and magnetic 
field fluctuations are correlated or not. In the case of negative <jR^ this question may easily be resolved: 
if the noises are correlated the linewidth is proportional to (Rd + crRd) 2 and may be significantly reduced 
if \Rd\ ~ |<ri?^|; if the noises are uncorrelated, the linewidth is proportional to i? 2 , + (aR^) 2 and the sign 
of uRd w iH not affect the linewidth. It could be interesting to study the possibility to achieve stable 
generation regimes for aRd < for practical FFOs, but as preliminary results demonstrate the curves of 
f2j(r) are similar to that presented in Fig. 5 for 7/0 = 0.5, i.e. they have positive derivatives. However, 
the optimization of practical FFO designs with respect to working parameters is out of scope of the paper 
and will be presented elsewhere. 
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V. CONCLUSIONS 



In the present paper the approach for description of the dynamics of dense soliton chain in long 
Josephson junctions for the case of large bias current and magnetic field is outlined. The approach allows 
to perform the detailed analytical study of flux flow oscillators with spatially inhomogeneous driving and 
to find an optimal profile of the bias current distribution in order to maximize the output power and to 
minimize the linewidth. Practically important example of FFO with unbiased tail has been considered 
and the existence of an optimal length of the unbiased tail has been demonstrated. The possibility to 
significantly reduce the linewidth by the design of FFOs with anticorrelated bias current and magnetic 
field fluctuations has been discussed. 
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FIG. 1: Current-voltage characteristic of FFO with homogeneous driving. Numerical solution of the sine-Gordon 
equation is presented by crosses, the fourth order approximation is given by solid line and the second order 
approximation is given by dashed line for a = 0.1, L = 5, T = 3. Inset: Enlargement of the IVC for flj < 1.2, the 
step (around flj « 0.95) missed in the second order is reproduced in the fourth order approximation (dimensionless 
units). 
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FIG. 2: Current-voltage characteristic of FFO with homogeneous driving for a = 0.04, L = 5, F — 5; crosses - 
computer simulations, solid line - theory (dimensionless units). 
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FIG. 3: Current-voltage characteristic of FFO with unbiased tail, a — 0.3, L — 5; for xo = 0, T = 5 - solid line; 
for xo = 2.5, T = 5 - triangles; for xo = 2.5, r = —5 - circles (see explanation in the text). Inset: IVC for a = 0.3, 
L = 5, xo = 2.5, T = —5; solid line - computer simulations, crosses - second order approximation (dimensionless 
units). 
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FIG. 4: Current-voltage characteristic of FFO with unbiased tail, a = 0.3, L = 5, V = —5, bias is from xo to L; 
solid lines - computer simulation, dots - theory. It is seen, that there is an optimal length of the bias tail both 
from the point of view of Rd and power of the main harmonic. Inset: IVC minus ohmic line for a — 0.3, L — 5, 
T = — 5 and xo = - open circles, xq — 2.5 - filled diamonds and xq — 4 - open triangles (dimensionless units). 
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FIG. 5: Voltage as function of magnetic field F for a — 0.3, L = 5 and different values of bias current: 1 - rjo = 0.5, 
2 - r)o = 1, 3 - r)o = 1.5; open circles - xo = 0, filled diamonds - xo = 2.5, open triangles - xo = 4; crosses - theory 
for r]o = 1, xo = 2.5 (dimensionless units). 
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